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We consider two-band double exchange model and calculate the critical temperature in ferromag- 
netic regime (Curie temperature). The localized spins are represented in terms of the Schwinger- 
bosons, and two spin-singlet Fermion operators are introduced. In terms of the new Fermi fields the 
on-site Hund's interactions are in a diagonal form and one accounts for them exactly. Integrating 
out the spin-singlet fermions we derive an effective Heisenberg type model for a vector which de- 
scribes the local orientations of the total magnetization. The transversal fluctuations of the vector 
are the true magnons in the theory, which is a base for Curie temperature calculation. The critical 
temperature is calculated employing the Schwinger-bosons mean-field theory. While approximate, 
this technic of calculation captures the essentials of the magnon fluctuations in the theory, and for 
2D systems one obtains zero Curie temperature, in accordance with Mermin- Wagner theorem. 
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PACS numbers: 75.47.Lx, 71.27.+a, 75.30.Ds, 75.50.Cc 



I. INTRODUCTION 

The manganites are ones of the prominent represen- 
tatives of strongly correlated systems, where the effect 
of correlations among electrons plays a crucial role. The 
growing interest on manganites is mainly due to observa- 
tion of resistivity changes by many orders of magnitude 
upon the application of small magnetic fields, an effect 
that carries the name of " Colossal Magnetoresistance" . 

The attempts to properly describe the physics of these 
materials by means of simple pcrturbativc methods typ- 
ically fail [l|]. One have to develop technic of calculation 
which capture the essential features of the compounds. 
The double exchange model is a widely used model for 
manganites [l|, [3, Q. In isolation, the ions of Mn have 
an active Sd-shell with five degenerate levels. The de- 
generacy is presented due to rotational invariance within 
angular momentum I = 2 subspacc. The crystal environ- 
ment results in a particular splitting of the five c?-orbitals 
{crystal field spliting) into two groups: the Cg and t2g 
states. The electrons from the Cg sector, which form a 
doublet, are removed upon hole doping. The t2g elec- 
trons, which form a triplet, are not affected by doping, 
and their population remains constant. The Hund rule 
enforces alignment of the three t2g spins into a 5 = 3/2 
state. Then, the t2g sector can be replaced by a localized 
spin at each manganese ion, reducing the complexity of 
the original five orbital model. The only important in- 
teraction between the two sectors is the Hund coupling 
between localized t2g spins and mobile Cg electrons. A 
strong impact on the physics of manganites has the static 
Jahn- Teller distortion which leads to a splitting of the de- 
generate Cg levels. 

The double exchange model has a rich phase diagram, 
exhibiting a variety of phases, with unusual ordering in 
the ground states. The procedures followed to obtain 
the phase diagram in one band model are different: nu- 
merical studies [H , dynamical mean field theory Q , and 
analytical calculations 0] , but four phases have been 



systematically observed: (i) antiferromagnetism (AF) at 
a density of mobile electrons n = 1, (ii) ferromagnetism 
(FM) at intermediate electronic densities, (iii) phase sep- 
aration (PS) between FM and AF phases, and (iv) spin 
incommensurable (IC) phase at large enough Hund cou- 
pling. The competition between spin spiral incommen- 
surate order or phase separation and canted ferromag- 
netism is also a topic of intensive study [1, 0i 9 • The 
phase diagram becomes more rich if the orbital degener- 
acy is accounted for [§, [l^l ■ 

Double exchange model is also successful in calculat- 
ing the critical temperature in the ferromagnetic regime. 
Predictions about Curie temperature in one band dou- 
ble exchange model are made using Monte Carlo tech- 
nique [1], High-temperature series expansion [Tl| . Dy- 
namical Mean-Field Approximation [131 and standard 
Mean-Field approach [ij]. The most striking feature of 
the critical temperature as function of fermion density is 
the symmetry with respect to n = 0.5. The Curie tem- 
perature is maximized at that point, and maximal value 
is different within different approaches. 

Very recently the critical temperature was calculated 
using two-band model with account for the Jahn- Teller 
effect [13, Again the characteristic feature of the 

temperature curves as a function of charge carrier density 
is the symmetry with the respect to n = 1. The new 
assertion is that the Curie temperature increases with 
increasing the interband hopping 
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It is impossible to require the theoretically calculated 
Curie temperature to be in accordance with experimen- 
tal results. The models are idealized, and they do not 
consider many important effects: phonon modes, several 
types of disorder. Coulomb interaction, etc. Because of 
that it is important to formulate theoretical criteria for 
adequacy of the method of calculation. In our opinion 
the calculations should be in accordance with Mermin- 
Wagner theorem p^ . It claims that in two dimension 
there is not spontaneous magnetization at non-zero tem- 
perature. Hence, the critical temperature should be equal 
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to zero. It is well known that Monte Carlo method of 
calculation does not satisfy this criteria [l| . It is difficult 
within Dynamical Mean-Field theory to make a difference 
between two dimensional and three dimensional systems. 
DEM is a good approximation when the dimensionality 
goes to infinity. This made the analytical methods im- 
portant even for the assessment of the numerical results. 

The paper is organized as follows: In Sec. II, start- 
ing from two band double exchange model, we derive 
an effective Heisenberg like model in terms of vector de- 
scribing the local orientations of the total magnetization. 
The transversal fluctuations of the vector are the true 
magnons in the theory. This is a base for Curie tempera- 
ture calculation. Sec. Ill is devoted to phase diagrams of 
the model in space of Hund's constant and charge carrier 
density. We calculate spin-stiffness constant as a function 
of density which is important step towards understanding 
the Curie temperature behavior as a function of charge 
carrier density. The results for the Curie temperature are 
reported in Sec. IV. A summary in Sec. V concludes the 
paper. Spin-stiffness constant calculations are presented 
in the Appendix. 



with an action given by the expression 
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where P is the inverse temperature and the Hamiltonian 
is obtained from equations (jll.ip and (jll.2p replacing the 
operators with the functions. 

We introduce spin-singlet Fermi fields 
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[(pil{T)Cil2{T) - ipi2{T)cill{T)] (II. 6) 



which are U{1) gauge variant with charge -1 and 1 re- 
spectively 

vf'f (r) = e-»^(^)vl/^(r) M/'f (r) = e-^(^)vl,^(r) (II.7) 



The equations pi.5|) and (|II.6p can be regarded as a 
SU(2) transformation: 



II. EFFECTIVE MODEL 
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We consider a two-band double-exchange model, with 
Hamiltonian: 



with = ^iii and = ^ii2. 

In terms of the new Fermi fields, electron creation and 
destruction operators have the form: 
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where 1,1' are band indexes, ij are site labels, a are the 
spin indices, cj^. and cua are creation and destruction 
operators for mobile electrons, /i is the chemical poten- 
tial. Sii is the spin of the conduction electrons, and 
is the spin of the localized electrons. The sums are over 
all sites of a three-dimensional cubic lattice, and 
denotes the sum over the nearest neighbors. We denote 

^11 = il, ^22 = ^2, and ti2 = ^21 = t' . 

In terms of Schwinger-bosons {ipiaT^p\„) the core spin 
operators have the following representation: 
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and the spin of the conduction electrons Sa is 

<i = l4.r^.'^a.' = {or^tiXa'^^^i^', (11.10) 



where 



(11.11) 



It is convenient to introduce three basic vectors which 
depend on the Schwinger-bosons 
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with the Pauli matrices (r^ ^t^ ,t^). 

The partition function can be written as a path integral 
over the complex functions of the Matsubara time Lpicy (r) 
((^^(r)) and Grassmann functions cuair) (c^^. ('''))■ 
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where Tf = ^S^. Then, the spin of the electrons can be 
represented as a linear combination of three vectors Sj, 
P, =T]+ iT? and P+ = - iT? 
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The basic vectors satisfy the relations = s^, Pf = 
P+2 = Si • Pj = S j ■ P+ = 0, and P,+ • P, = 2. Using the 
expression ()IL13P for the spin of itinerant electrons, the 
total spin of the system 



- S, + s,i 
can be written in the form 



Si2 



(11.14) 
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^f+^Ulp+J2^t+^f^ (11.15) 



The gauge invariance imposes the conditions 

) = (*,^+*^) = 0. As a resuh, the dim 
less magnetization per lattice site ((5*'°*)^) reads 



(11.16) 

Let us average the total spin of the system (Eq. III.15[) 
in the subspace of the itinerant electrons (S*°')/ = M^. 
The vector identifies the local orientation of the total 
magnetization. Accounting for the gauge invariance, one 
obtains the following expression for Mi 



where 
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can be thought of as an "effective spin" of the system 
(M? = M^). Now, if we use Holstcin-Primakoff repre- 
sentation for the vectors Mo 
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MJ = Mji - iMj2 = aj ^2M ~ aja, (11.19) 
A'/f = M - a+a, 

J J 

the bose fields aj and are the true magnons in the 
system. In terms of the true magnons the Schwinger- 
bosons pi.2p have the following representation 
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Replacing in Eqs. pi.sp and pi.6[) for the spin-singlet 
Fermions and keeping only the first two terms in \/M 

expansion \/l — jIT '^t'^i — 1 ~ tm a^ai + . . . we obtain 
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The equations pi.2ip and pi.22|) show that the singlet 
fermions are electrons dressed by a virtual cloud of re- 
peatedly emitted and reabsorbed magnons. 

An important advantage of working with A and B 
fermions is the fact that in terms of these spin-singlet 
fields the spin-fermion interaction is in a diagonal form, 
the spin variables (magnons) are removed, and one ac- 
counts for it exactly: 
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Replacing all this into the action (|TL4|, we can rewrite 
it as a function of the Schwinger-bosons and spin-singlet 
fermions. The resulting action is quadratic with respect 
to the spin-singlet fermions and one can integrate them 
out. The effective Hamiltonian of the theory, in Gaussian 
approximation, is given by: 

/leff = P E i^-to-^ + flj^^j ^ O't'^'j ^ aljUi) (11.24) 



where p is the spin stiffness (|A.42p . Detailed calculation 
are given in the appendix. Based on the rotational sym- 
metry, one can supplement the Hamiltonian (Eq. III.24[) 
up to an effective Heisenberg like Hamiltonian, written 
in terms of the vectors Mi 



(H.25) 



where J = p/M . The ferromagnetic phase is stable if the 
effective exchange coupling constant is positive J > 0. 

It is important to highlight the difference between 
the above effective theory pi.25p and Ruderman-Kittel- 
Kasuya-Yosida (RKKY) theory. The RKKY effective 
Hamiltonian is an effective Heisenberg like Hamiltonian 
in terms of core spins Si, obtained averaging in the sub- 
space of the itinerant electrons [l3|- The subtle point 
is that if we use a Holstein-Primakoff representation for 
the localized spins Si , the creation and annihilation bose 
operators do not describe the true magnon of the sys- 
tem [Tsj . The true magnons are transversal fluctuations 
corresponding to the total magnetization which includes 
both the spins of localized and delocalized electrons. 
Therefore the RKKY validity condition requires small 
Hund's coupling, and small density of charge carriers, 
which in turn means that the magnetization of the mo- 
bile electrons is inessential. In contrast of RKKY theory 
the effective model pi.25p is written in term of vectors 
Mi which describe the local orientations of the total mag- 
netization, and the bose operators in pi.l9p are the true 
magnons in the theory. This is essential when one calcu- 
lates the Curie temperature. The effective model pl.25p 
is obtained integrating out the spin-singlet fermions pI.SP 
and pi.6p . In terms of these fermions the Hund's inter- 
action is in a diagonal form and we account for it exactly. 
Hence, the effective theory pi.25p is valid for arbitrary 
values of Hund's constants and for all densities of charge 
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III. PHASE DIAGRAMS 

Here we illustrate some of the features of our model, 
namely the phase diagrams and how they change when 
we vary the model's parameters. Since calculating is 
closely related to calculating spin stiffness p, we have ex- 
amined the behavior of p as a function of electron density 
in details. 

The physics of the model depends on dimcnsionlcss pa- 
rameters Ji/ti, J2/ii, ti and t' /ti. Throughout this 
chapter we fix the scale setting = ^2 = 1, and use Ji, 
J2 and t' as a free parameters of the model. Also in this 
section t' — 0.5, and since we are describing manganite 
materials, we have set s = 3/2. The phase diagrams on 
FigH] are constructed by plotting the curve p = in co- 
ordinate system of carrier density n and J2/ti for fixed 
ratio Ji/J2- Regions where p > correspond to the 
ferromagnetic phase (FM), while those with p < Q are 
denoted here simply as non-FM ones, since describing all 
possible phases is not the purpose of this paper. 
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FIG. 1: Phase diagrams for the four cases under consid- 
eration: a) Ji = J2 ~ 15, b) Ji = 16.5, J2 — 15, c) 
Ji = 17.5, J2 = 15 and d) Ji = 25, J2 = 15. White regions 
correspond to FM phase, grayed ones to FM instability. 




FIG. 2: (color online) p as a function of n for the forth case 
J1/J2 = 25/15, corresponding to fig. [TJi. 



0.08 



0.04 



0.02 



1 1 1 1 1 1 1 


i i 1 i i i i 1 i 1 i 1 

J,/J2 = 17.5/15 ■ 




J, = 45 








J^ = 30 




\'-, \ J, = 25 . 




\ ~\ \ ^2 = 20 - 




\ \ \ -^2=^^ ■ 


1 . 





0.5 1.0 1.5 2.0 

n 



FIG. 3: (color online) p as a function of n for the third case 
J1/J2 — 17.5/15, corresponding to fig. [!};. 



We consider four different cases for the ratio J1/J2, 
namely J1/J2 = 15/15, where the bands are degenerated, 
and three cases with increasing bands' splitting J1/J2 = 
16.5/15, J1/J2 = 17.5/15 and J1/J2 = 25/15. 

When the bands are split, we observe an island of fer- 
romagnetic instability around the line n = I. The lowest 
point of the island is denoted by . Increasing the ratio 
J1/J2 increases the island by lowering the value of 
(see FiglT]). For the degenerated bands (J1/J2 = 1), we 
have not observed the island of ferromagnetic instability, 
up to values as large as J2 = 300. 

The value of the spin stiffness constant p depends on 



the point (71, J2) in the phase diagram, for fixed ratio 
J1/J2. When the point approaches the boundary of the 
ferromagnetic phase, the value of p decreases and reaches 
zero on the boundary (by the definition of the boundary) . 
We have calculated the spin stiffness constant p as a func- 
tion of n, for fixed J2 at zero temperature. There are 
three distinctive cases: J2 > J2 7 J2 ^ J2 and J2 ^ J2. 

For the first one, the presence of ferromagnetic insta- 
bility near n = 1 results in a function p{n), which consists 
of two pieces, one for n in the interval (0, 1), and another 
for n in the interval (1,2). For the ratio J1/J2 = 25/15 
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(Figl2]), such are the curves corresponding to J2 = 22, 
J2 = 15 and J2 = 10; for J1/J2 = 17.5/15 (FigO the 
curve corresponding to J2 ~ 45; and for J1/J2 = 16.5/15 
(FigH]) the one corresponding to Ji = 75. The impor- 
tant characteristic of aU these curves is the existence of 
two maxima, one within interval (0, 1), and another one 
within interval (1, 2). The global maximum is within the 
interval (0, 1), which is result of the presence of two phase 
boundaries in the other interval. In the case of degener- 
ated bands, the absence of island of instability leads to 
the absence of such type of function pin) . 

For the second one, J2 is very close to J|, hence near 
n = 1 the spin stiffness constant is very small. As a result 
p{n) is a function with two maxima and one minimum. 
For the ratio J1/J2 = 25/15 (FiglJ), such is the curve 
corresponding to J2 = 7; for J1/J2 = 17.5/15 (Figl3]) 
the curves corresponding to J2 = 30 and J 2 = 25; and 
for J\l Ji = 16.5/15 (Figdl) the ones corresponding to 
J2 — 60 and J2 = 50. The minimal value of the function 
decreases when J2 approaches . Again, in the special 
case of degenerated bands, there is no such a curve. 
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FIG. 5: (color online) p as a function of n for the degenerated 
case J1/J2 = 15/15, corresponding to fig. [1^. 
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FIG. 4: (color online) p as a function of n for the second case 
J1/J2 = 16.5/15, corresponding to fig. [TJj. 

For the third case, J2 <C J|, we are well below the 
island of ferromagnetic instability, and the function p{'n) 
has only one maximum within the interval (0,1). For 
the ratio J1/J2 = 25/15 (FiglJ), such is the curve cor- 
responding to J2 = 5; for J1/J2 = 17.5/15 (Figl3]) the 
curves corresponding to J2 = 20 and J2 = 15; and 
for J1/J2 = 16.5/15 (FigHI) the ones corresponding to 
J2 — 30 and J2 = 15. In the case of degenerated bands, 
for all values of J2 the curves are of this type, because of 
the absence of instability island (Fig[5]). 

It is widely known fact, that near n = 2 the system is 
ferromagnetically unstable and the spin stiffness constant 
approaches zero when the carrier density approaches two. 



One can see from the phase diagrams (Fig[T]), that with 
decreasing of J2 , the ferromagnetic instability sets in for 
smaller value of n. As a result wc obtain that with de- 
creasing J2, the value of n for which p[n) — decreases. 
This is best observed in the case of degenerated bands, 
where we have examined broader set of values for J2. 
With increasing of J2, the point p = gets closer to 
n = 2 and the curve becomes more symmetric. 

Since the maximum value of the Curie temperature 
corresponds to the maximum value of the spin stiffness 
constant, it is important to see how this value depends on 
the ratio J1/J2 for fixed J2. We choose the case J2 = 15 
(red lines in the figures), to compare our results with the 
results in p^ . Increasing J1/J2 increases the ferromag- 
netic instability, see phase portraits Fig. [l] lowering the 
point J| which in turn leads to decrease of the maxi- 
mum value of p and qualitative changes in its behavior 
as a function of density n. 



IV. CURIE TEMPERATURES 

To calculate the Curie temperature we utilize the 
Schwinger-bosons mean- field theory [l^, . The advan- 
tage of this method of calculation is that for 2D systems 
one obtains zero Curie temperature, in accordance with 
Mermin- Wagner theorem [l6|. So, while approximate, 
this technique of calculation captures the essentials of 
the magnon fluctuations in the theory. 

To proceed, we represent the effective spin vector 
by means of Schwinger bosons (f>ia , 4>f^ , 

Mr = ^ E <^-' '^-'^^ = 2M (IV. 1) 
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Next we use the identity 



and rewrite the effective Hamiltonian in the form 
J 



(IV.2) 



<i]> 



where the second term in pV.2[) is equal to the constant 
M^, because of the constraint (|IV.1[) . and we drop it. To 
ensure the constraint (jIV.ip we introduce a parameter (A) 
and add a new term to the effective Hamihonian (jfV.3[l 



(IV.4) 



We treat the four-boson interaction within Hartree- 
Fock approximation. The Hartree-Fock hamiltonian 
which corresponds to the effective hamiltonian (jIV.4p 
reads 



J 



<»i> 



E 



E[ 

<lj> 

Xj2{^tAa-2M) (IV.5) 



where (uy ) arc Hartrec-Fock parameters to be de- 
termined self-consistcntly. We are interested in real 
parameters which do not depend on the lattice sites, 
Uij ~ u. Then in momentum space representa- 



tion, the Hamilonian pV.5|) has the form 

hn-F = Y^"' " ^^^^^ + E ^'^'f't^k (IV.6) 

k 

where N is the number of lattice sites and 

Ek ^ X — Ju (cos kx + cos ky + cos fc^) (IV. 7) 

is the dispersion of the 0fe-boson (spinon). 

The free energy of a theory with Hamiltonian pV.6p is 
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where T is the temperature. The self-consistent equa- 
tions for parameters u and A are 



We obtain a system of two equations 
2 1 



(IV.9) 



-— ^nfc {coskx + cos ky +cosk^) (IV. 10) 
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where is the bose function 
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(IV.12) 



To ensure correct definition of the bose theory (|IV.6|) 
we have to make some assumptions for the parameter 
A. For that purpose its convenient to represent it in the 
form 



A = 3 Ju + fiJu 



(IV.13) 



In terms of the new parameter fi, the bose dispersion is 



Ju (3 — cos kx — cos ky — cos k^ + fi) (IV. 14) 



and the theory is well defined for fi > 0. 

We find the parameters fj, and u solving the equa- 
tions pV.lOmV.lip . For high enough temperatures both 
n{T) > and u{T) > and the excitation is gapped. It 
is the spinon excitation in the theory in the paramagntic 
phase. Decreasing the temperature leads to decrease of 
/i(r). At temperature Tq it becomes equal to zero and 
long-range excitation emerges in the spectrum. Hence 
the temperature at witch /i reaches zero is the Curie 
temperature. We set /i = in the system of equations 
PV.10IIIV.11[) and obtain a system of two equations for 
the Curie temperature Tc and the parameter u, which 
is the renormalization of the exchange constant at Curie 
temperature. 
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cos kx + cos ky + cos kz 



3 iV ^ — ^ — (3— COS A:^ — COS fc^ — COS A:^ ) -, 



M 
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(IV.15) 



A'^ ^ ^ TTW' — (3 — COS fc^ — COS /i:^ — COS fc^ ) -, 



To calculate Tc we solve the above system of equations, 
with M(r) and p(r) calculated from equations pi.l8|) 
and (|A.42[) respectively. 

The results for the Curie temperature Tq as a function 
of charge density n are plotted on figures [Sjni for two dif- 
ferent values of interband hopping parameter t' . In all 
cases we have set ti = <2 = 1 , and consider a theory with 
s = 3/2 for the core spins. We want to compare our re- 
sults with the results in [l^, so we have set Ji = 15 and 
t' = 0.5. Since our method of calculating Tc involves 
p{T), the resulting curves are very similar to the ones in 
section HI (Figl2]15]), and bear their characteristics. As 
above we have three different choices for the parameters: 
J2 > J2 ^ Ji and J2 ^ ^2- For the biggest ratio we 
consider J1/J2 = 25/15 (see Figl6|), we have a two-piece 
function, since our chosen value J2 > J2 ■ With decreas- 
ing the ratio, the point J| moves to higher values and J2 
becomes smaller than J|. As a result the curves we ob- 
tained have only one maximum. It is important to note, 
that even for the degenerated case J1/J2 = 1, the point 
at which Tc reaches zero is smaller than n ~ 2 unlike 
the results in [l3|. This can be also seen from the phase 
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FIG. 6: (color online) Curie temperature Tc as a function of FIG. 8: (color online) Curie temperature Tc as a function of 
charge density n for Ji = 25, J2 = 15, with t' = 1/2 and charge density n for Ji = 16.5, J2 = 15, with t' = 1/2 and 
t' = l t' ^1 
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FIG. 7: (color online) Curie temperature Tc as a function of 
charge density 7i for Ji = 17.5, J2 — 15, with t' — 1/2 and 
t' = l 



FIG. 9: (color online) Curie temperature Tc as a function of 
charge density n for the degenerated case Ji — 15, J2 = 15, 
with t' = 1/2 and t' ^ 1 



portraits (Fig[T]), where the instability of ferromagnetism 
near the charge carriers density n = 2 is evident. One 
can also note that decreasing the ratio J1/J2 leads to 
increase in the maximum value of Tc{n). 

The most important difference between our curves of 
critical temperature as a function of n and the ones in [l3 | 
is the lack of symmetry with respect to n = 1. This asym- 
metry originates from asymmetry in the phase diagrams, 
which in turn leads to asymmetry of the spin-stiffness 
curves. 



We have also examined the effect of t' . Increasing its 
value results in enlargement of the ferromagnetic insta- 
bility island around n = 1. The value of J| decreases 
and the width of the island increases. This in turn leads 
to both quantitative and qualitative changes in Tc(ri) 
(see the black lines of Figl6][9]) . An important conclusion 
is that Curie temperature increases with increasing t' if 
the bands are strongly split (Figl5][7]), while for weakly 
split bands we obtained an opposite behavior, the critical 
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temperature decreases (Figl8]|9l). 



APPENDIX A: CALCULATION OF p 



V. SUMMARY 

In summary, we have calculated the Curie tempera- 
ture in two-band double exchange model. First we re- 
duced the model to an effective Heisenberg type model 
for a vector which describes the local orientations of 
the total magnetization. Next, we use Schwingcr- 
bosons mean-field theory to calculate the critical temper- 
ature. This technique of calculation is in agreement with 
Mermin- Wagner theorem, which means that employing 
our method of calculations for 2D system one obtains 

There are many quantitative and qualitative differ- 
ences between our results and the results obtained within 
Dynamical Mean-Field Theory or Monte Carlo simula- 
tion approach. Maybe the most significant difference is 
that DMF and MC calculations lead to temperature's 
curves, as a function of fermion density, symmetric with 
respect to n = 1. This is not the result in the present 
paper. The asymmetry of the Tc{n) curves with respect 
to n = 1 is a characteristic feature in our approach. This 
asymmetry is seen looking at the phase diagrams. It 
predetermines the asymmetry of the spin-stiffness curves 
which lead directly to the asymmetric Tc{n) curves. The 
symmetry mentioned in the paper is possible if the 
two bands are degenerated and Hund's constants are un- 
physically large. 

Another important result reported in previous papers 
is that the critical temperature increases with increasing 
the interband hopping. Our calculations show that this is 
true when the band splitting is strong. If the bands split 
weakly, the assertion is opposite. The Curie temperature 
increases when the interband hopping decreases. 



Here we present a detailed derivation of the spin stiff- 
ness constant p. We start from the two-band hamiltonian 
(|IlT|) and rewrite it in terms of Scwinger bosons (|II.2[) 
and spin-singlet fermions (jll. 51111.6]) . 

The resulting action is quadratic with respect to the 
spin-singlet fermions and one can integrate them out. To 
do so, it is convenient to represent the action as a sum 
of three terms: 



where Sf is the free fermion action: 



(A.l) 



S,=Jdrl^(^+^-^^,,, + ^+^^^,2)j+Hj (A.2) 



with free fermion hamiltonian Hf 



il ila 



<ij><T 



(A.3) 



<ij><y 



<ij>(T 
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We remind of the notations '^f m = and 'if ii2 = , so 
that the sum over a in the above equation is a sum over A 
and B. It is convenient to represent the term describing 
spin-fermion interaction as a sum of two terms: 



J 



Si% = Jdr 





2.S 



2s 



<i]> 



2s) (*i^*f, +*+^*fJ 



2s 



<i]> 



2i) 



(A.4) 



'2s 



<ij> 



(A.5) 



To diagonalizc the free fcrmion Hamiltonian Hi wc 
switch to momentum space, 



Ik 

+ E (*^c.*2fe. + (A.6) 

where the dispersions for A and B fermions are 
e^i = -2ti ^ cos fc^ - sJi- fi 

k 

^kl = -2t; ^ COS fc^ + S J; - fi ^^ j-^ 

k 

e(k) = -2t'^cosfc^ 



The Hamiltonian is diagonahzed by means of the trans- 
formation: 



Solving the equations for u and v gives 



(A.8) 



,R 



z;|^ = sign(e(fc))\/^(l-x«) 



with ^k — — 7 

The resulting Hamiltonian is 



R _ R 

^2k ^Ife 



'2k 



(A.9) 



(A.IO) 



-e?ky 



^IkJlk Jlk + ^Ife/lfc /ife 



B J-+-B j-B 



pA r+A fA , pB f+B fB 
^2kJ2k J 2k ' ^2kJ2k J 2k 



with dispersions for the quasi-particles 



E 



Ik 



E2k = 



2 



(A.12) 



^B ^ ^2k + ^ffc 

1'' 2 2 



rB _£B^ 



-B I ^B 



(A.13) 



1 



To write the spin-fermion interaction Sg-t in terms of 
the new fermions we introduce the notations 



f+ — ( f +-4 f+B f+A f+B \ f 
Jki ^ \Jlki Jlki J2ki J2ki ) J ^2 



Now we rewrite the action in the form 



fA 
J2k2 

fB 
yJ2k2 J 



(A.14) 



Ss-{ - J dndr^Y^ft {ri)Wk,k2 (T1-T2) /fe,(T2) (A.15) 
n ^1 ^2 



where the matrix Wk-^k2 {'^1 — ^2) is a sum of two terms 



Wk,k2 (ri-T2) = Wl^^_ {T^-T2)+WtX, {T1-T2) . (A.16) 
(A. 11) ^kik2 ('''1 ~ ''2) is the free fermion action 
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iSr + E: 



dT2 ^ lk2 





V 



^ + ^2^ 



\ 






while W^fcjfcj (ti — T2) is a sum of two matrixes , corresponding to S\_i^ and S\_!_^ 



( Ku,k2 (t2) 
Lk,k2 {t2) 
Nk,k2 {t2) 



\ 






V 







Ok,k2 (t2) / 



+ 5(ti-T2) 



/ Afc^fe. (T2) Bfc^fc, (t2)\ 

-Bfc.fc, (T2) -Ffcifc. (T2) 

VGfe.fc, (t2) /fcifc, (t2) / 



with matrix elements 
Kk^k2 (7-2) = 

ifcifea ('^2) = 

^felfc2 (T2) = 



2sN ^ \ ^ 3 

fc \fj.=l / <ii> 



2s TV ^ \ ^ 3 

fc \^=l / <ij> 



2.S iV ^ \ ^ 3 



2s A/' ^ \ ^ 3 

fe \M=1 / <ij> 



(A.17) 



(A.18) 



(A.19) 



Ak.k2 (r2) = E i^M2 - ^ti^t2) (e-*'^^+'^-- - + - " 

(A.20) 

171 / \ 11 / + + + +^ / -ikiri+ikon -ikiri+iknri\ V, A B , A B ./ABi.rABl 

Ek,k2 {t2) = -—^^2^ [fTi'Pj2 - fji^72) (e ^ ^+ ^ • - e 1 ^ ^) [iiWfc^Ufc^ - t2U^,v^^ - t v^.v^, + t u^^u^^l 

<i]> 



-ikirj-\-ik2ri —ikiri-\-ik2rj 



) + i2« + + t'« 



<ij> 
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<ij> 

<ij> 



(A.21) 



-ikiVi -}-ik2rj 



r 



where is the number of lattice's sites. Integrating the 
fermions out we obtain the effective action Scs 

S^ff = - In det = -Tr hi W ( A.22) 

Using the properties of the logarithm 

Tr In = Tr In (1^° + W"'^) = Tr In W° 

+Trln (l + {W°y^ M^'"*) (A.23) 



where the term TrlnW^" doesn't depend on the 
Schwinger bosons and we have dropped it. Finally, we 
expand the effective action in powers of 

Xk,k, {ruT2)^Y. jd^Wl, (Ti,r)]-V^f,(r,r2) (A.25) 



The result is 



we rewrite the effective action in the form 



S^s = -TvX + -TrA^ 
2 



(A.26) 



S-eff^-Trln 1 + (VK°) W 



(A. 24) The inverse matrix (W^^^/^j) ^ is given by 



( Kk.Stu,{Tl-T2) 

4,fe.^2l(ri-T2) 

V 



(A.27) 



^k^k2S2u, (n -7-2) / 



where S','^^ {u) 



E, 



(a ^ A OT B, I = I or 2). 



Ik 



Replacing (jA.lSp into (|A.25p . we end up with two terms 
for Xfejfej, one which is diagonal XP'l , and one with 



(2) 

zero diagonal elements A), , . Hence, one obtains for the where 



k\k 



trace of the matrix X 

Tr A = TrX(i). 



(A.28) 



A 



(1) 



(Ti,r2): 



/5il(ri-r2)iffc,fc,(r2) 

5il(ri-r2)ifc,fe,(T2) 


V 



qA 



\ 



{ri-T2)Nk,kAr2) 

Sg^{n~T2)Ok,kAr2)J 



(A.29) 
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and the first term in the effective action (jA.26p is 



'eft 



+tl 



1 1 

'2^N 



k \m— 1 







t2 



{nt) n^k + {vtYnt + {uiY n^, + (.f ) nf. 



■ (A.30) 



To calculate the contribution of 5^^^ to the spin- and keep the terms quadratic with respect to the 
stiffness constant p in (jll.24p we use the Holstcin- magnons , . The result is 
Primakoff representation for the Schwinger bosons 



, * + / * + 



(A.31) 



Pi = 



1 1 



2M N ^ , 

A: \m=1 



ECOS k,. 



2t' 



) + ^kVk « - O 



■2k ~ "Ifc 



+tl 



{utYnt+{v^Y nt+i4)"nU{viY < 



+ t2 



{utYnt+{v^Y<+i^kYr^2k + {viY< 



(A.32) 



Calculating the contribution of the second term in magnons, hence it doesn't contribute. The contribution 
(|A.26| to the effective hamiltonian (|II.24p we account comes from S'^'^'f ((XS)) which leads to the matrix AT'^). 
for the fact that matrix is quadratic with respect to 
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/ ^il(Tl-T2)Afe,fc,(r2) ^il(Ti-T2)i?,,fe,(r2)\ 

^iKri-T2)Cfe,fe,(r2) S^,^in-T2)Dk,kAr2) 

S^,An^T2) Ek,kAr2) ^2l(n-r2)F;.,;.,(r2) 

V^2KTi-r2)Gfe,fe,(r2) 52l(n-r2)4,.,(r2) J 

I 



After some algebra we arrive at the following representa- tion of the second term in (jA.26p 
I 



^iff = / dndT2Y, [sX (^1 - ^2) Ak,k, {t2) 5f,^ {t2 - ri) C^^k, (n) + S^, (n - T2) B,,,, (rs) ^2! (^2 - n) G^^fc, (n) 

(A.34) 



felfc2 



+^2l-, (n - r2) Ek,k, {t2) 5ffe, (r2 - ri) Dfe.fc, (n) + 52I,, (n - r2) F^.^, (r2) S^,.^ {t2 - n) 4,^, (n) 

I 



Switching from imaginary time r representation to fre- 



quency Lu representation we calculate the expressions in 
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small Lo approximation. The result is a 

(A.35) 

Next we make a change of wave-vectors variables fci = 



q+^k, ki = q—^k, and calculate the expressions in small 
wave-vector k approximation. The expression (|A.34p cal- 
culated in small frequency and small wave- vector approx- 
imation has the form 



'off 



duj 1 
2^N 



EE 




(A.36) 



,B 



'Uq + t'Uq 



Our third step is to express the products of the Green 
functions, in the above equation, in terms of the fermi 
hmction n{E) = l/(e^ + 1) 



2tt 



n {E^g) - n 



^2q - J^2q 



(A.40) 



duj 
2^ 



n (E^ 



1?/ 



(A.37) 



^2q ~ ^Iq 



Finally wc use the representation of the Schwinger bosons 
(|A.3ip . To calculate the contribution to the spin-stiffness 



/d n (E^ ) — n (E^ ) constant it is enough to keep only the quadratic terms 

^Siiuj)S^Juj) ^ —-^ (A.39) with respect to magnons 
Efq-E^^ 



P2 



2 1 
MV 



EE 



d . 2 ' 

sm 



+ {hvfuf ~ hufvf - t'vfvf + t'u^u^) 



2 ( n (SfJ - n [E 



WB zpA 



V E.fq - E^g 



2(n{E ?^)-n{E ^^) 

•B _ T^A 
Iq ^2q 



\ - E2 

2( n {Ei^)-n{E^^y 

\ E^q - E2q 



(A.41) 



The spin-stiffness constant in the effective action pi.24[) is a sum of the expressions (|A.32[) and (jA.41|) 
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ti 1 sr-^ I \r-^ COS fc^i 
2MV^ 1 2^ ^ 



\/^— 1 
/ d 



^2 1 / COS fcp 



k \l-L—l 



^•^^ / cos kf^ 



{r^2k - nt) + i-f [nl - «ffe) 



EE 



sm 



'„.A„.B I +i„.A„.B\ 



,A„,B I j-i„.A„.B I j-'„,A„ B\ 



(< 




\E?k 


^Ik 


f< 






- Efk 


(< 




W 




( < - 





~ E^ 



(A.42) 
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